Cosmological evolution of finite temperature Bose-Einstein Condensate dark matter 
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Once the temperature of a bosonic gas is smaller than the critical, density dependent, transi- 
tion temperature, a Bose - Einstein Condensation process can take place during the cosmological 
evolution of the Universe. Bose - Einstein Condensates are very strong candidates for dark mat- 
ter, since they can solve some major issues in observational astrophysics, like, for example, the 
galactic core/cusp problem. The presence of the dark matter condensates also drastically affects 
the cosmic history of the Universe. In the present paper we analyze the effects of the finite dark 
matter temperature on the cosmological evolution of the Bose-Einstein Condensate dark matter 
systems. We formulate the basic equations describing the finite temperature condensate, repre- 
senting a generalized Gross-Pitaevskii equation that takes into account the presence of the thermal 
cloud in thermodynamic equilibrium with the condensate. The temperature dependent equations 
of state of the thermal cloud and of the condensate are explicitly obtained in an analytical form. 
By assuming a flat Friedmann-Robertson- Walker (FRW) geometry, the cosmological evolution of 
the finite temperature dark matter filled Universe is considered in detail in the framework of a two 
interacting fluid dark matter model, describing the transition from the initial thermal cloud to the 
low temperature condensate state. The dynamics of the cosmological parameters during the finite 
temperature dominated phase of the dark matter evolution are investigated in detail, and it is shown 
that the presence of the thermal excitations leads to an overall increase in the expansion rate of the 
Universe. 

PACS numbers: 04.50.Kd, 04.20.Cv, 04.20.Fy 
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I. INTRODUCTION 



The Concordance Cosmological Model, usually re- 
ferred to as the ACDM (A+ Cold Dark Matter) model, 
has proved to be very successful in explaining cosmologi- 
cal observations across a wide range of length scales, from 
the cosmic microwave background (CMB) anisotropy to 
the Lyman-a forest \^. In this model, nonbaryonic col- 
lisionlcss cold dark matter makes up to 23% of the total 
mass content of the Universe. In the ACDM model dark 
matter consists of cold neutral weakly interacting massive 
particles, beyond those existing in the Standard Model 
of Particle Physics. However, up to now no dark mat- 
ter candidates have been detected in particle accelerators 
or in direct and indirect searches. Many particles have 
been proposed as possible candidates for dark matter, the 
most popular ones being the Weakly Interacting Massive 
Particles (WIMP) and the axions (for a review of the 
particle physics aspects of dark matter see [2[ ) . The in- 
teraction cross section of dark matter particles with nor- 
mal baryonic matter is assumed to be extremely small. 
However, it is expected to be non-zero, and therefore 
the direct experimental detection of dark matter parti- 
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cles may be possible by some existing or future detectors 
[J]. Superheavy particles, with mass > 10^° GeV, have 
also been proposed as dark matter candidates. But in 
this case observational results show that these particles 
must either interact weakly with normal matter, or they 
must have very heavy masses above lO""^^ GeV Q. Scalar 
field models, or other long range coherent fields coupled 
to gravity have been considered as possible candidates 
for galactic dark matter Q. The possibility that dark 
matter could be described by a fiuid with non-zero effec- 
tive pressure was also investigated [y, |7[. In particular, 
it was assumed that the equation of state of the dark 
matter halos is polytropic Q. The fit with a polytropic 
dark halo improves significantly the velocity dispersion 
profiles. The possibility that dark matter is a mixture 
of two non-interacting perfect fluids, with different four- 
velocities and thermodynamic parameters, was proposed 
recently in Q. It has been also suggested that galac- 
tic dynamics of massive test particles can be understood, 
without considering dark matter, in the context of mod- 
ified theories of gravity [ly . 

The experimental observation of the Bose-Einstein 
condensation of dilute alkali gases [ll[ represented a ma- 
jor advance in contemporary condensed matter physics. 
At very low temperatures, all particles in a dilute Bose 
gas condense to the same quantum ground state, form- 
ing a Bose-Einstein Condensate (BEG). In order for the 
condensation to occur particles must become correlated 



with each other so that their wavelengths overlap, that 
is, the thermal wavelength Ay must satisfy the general 
condition Xt > I, where I is the mean inter-particles 
distance. The condensation takes place at a tempera- 
ture T < 2TTh'^/mkBn^^^, where m is the mass of the 
particle in the condensate, n is the particle number den- 
sity, and ks is Boltzmann's constant [ij]- A coherent 
condensed state always develops if the particle density 
is high enough, or the temperature is sufficiently low. 
Quantum degenerate gases have been created in labora- 
tory by a combination of different laser and evaporative 
cooling techniques, opening several new lines of research, 
at the border of atomic, statistical and condensed matter 
physics [l3l • Recently, the Bose-Einstein condensation of 
photons has been observed in an optical microcavity jl3| . 

The possibility that dark matter, representing a sig- 
nificant amount of the total matter content of our Uni- 
verse, is in the form of a Bose-Einstein Condensate was 
analyzed in detail in [IJ]. By introducing the Madelung 
representation of the wave function, it follows that dark 
matter can be described as a non-relativistic, Newto- 
nian Bose-Einstein gravitational condensate gas, whose 
density and pressure are related by a barotropic equa- 
tion of state, which, for a condensate with quartic non- 
linearity, has the polytropic index one. The validity of 
the model was tested by fitting the predicted galactic 
rotation curves with a sample of rotation curves of low 
surface brightness and dwarf galaxies, respectively [ij- 
Il6l |. In all cases a very good agreement was found be- 
tween the theoretical rotation curves and the observa- 
tional data. Therefore dark matter halos can be de- 
scribed as an assembly of light individual bosons, oc- 
cupying the same ground energy state, and acquiring a 
repulsive interaction. This interaction prevents gravity 
from forming the central density cusps. The condensate 
particle is light enough to naturally form condensates 
of very small masses that later may coalesce. Different 
astrophysical properties of condensed dark matter halos 
have been extensively investigated [131 ■ 

The cosmological study of the Bose-Einstein conden- 
sate dark matter was initiated recently in [l^l and jl9| . 
respectively, with the study of the global cosmological 
evolution of the cold Bose-Einstein Condensate dark mat- 
ter, and with the analysis of the perturbations of the 
condensate dark matter. The obtained results show sig- 
nificant differences with respect to the pressureless dark 
matter model, considered in the framework of standard 
cosmology. Therefore the presence of condensate dark 
matter could have modified drastically the cosmological 
evolution of the early Universe, as well as the large scale 
structure formation process. The cosmological details of 
the Bose-Einstein Condensation process have been ana- 
lyzed for the case of the cold dark matter in [2^. The 
evolution of the cosmological inhomogcneities in the con- 
densed dark matter and the observational implications on 
the Cosmic Microwave Background (CMB) spectra have 
been considered in |2l| . 

All the mentioned previous researches have been done 



by assuming a zero temperature condensed dark matter. 
This assumptions is certainly a very good approxima- 
tion for the description of dark matter in thermodynamic 
equilibrium with the cosmic microwave background, and 
for the analysis of the galactic rotation curves. It is al- 
ready well established in condensed matter physics that 
the zcro-tcmpcrature Gross-Pitaevskii equation gives an 
excellent quantitative descriptions of the Bose-Einstein 
condensates for T < 0.5Ttr, where Ttr is the Bose- 
Einstein transition temperature [22|. This condition is 
obviously satisfied by the dark matter halos of the low- 
redshift galaxies. However, in the early Universe, im- 
mediately after the condensation, finite temperature ef- 
fects could have played an important role in the dark 
matter dynamics, and significantly influence the cosmo- 
logical evolution. The study of the finite temperature 
Bose-Einstein Condensate dark matter was initiated in 
[23j , where the first order temperature corrections to the 
density profile of the galactic halos and the static proper- 
ties of the condensates interacting with a thermal cloud 
have been obtained. 



It is the purpose of the present paper to study the cos- 
mological evolution of the finite temperature gravitation- 
ally self-bound Bose-Einstein dark matter condensates. 
As a first step in our study, by using the Hartree-Fock- 
Bogoliubov and Thomas-Fermi approximations, respec- 
tively, we obtain the equations of state of the arbitrary 
finite temperature dark matter condensate interacting 
with a thermal cloud. A finite temperature condensate 
can be described in terms of two fluids, the condensate 
proper, and the thermal excitations. Once the tempera- 
ture of the system decreases, the number of the thermal 
excitations drops, and at zero temperature the system 
consists of the condensate only. Therefore the dynamics 
of a finite temperature Bose - Einstein -Condensate can 
be described as a system of two interacting fluids. By 
using the specific equations of state of the thermal exci- 
tations and of the condensate, the cosmological dynamics 
of the Universe is considered within a two-fluid dark mat- 
ter model for a flat Friedmann-Robcrtson- Walker (FRW) 
type cosmological model. The physical parameters of the 
dark matter are obtained by numerically solving the evo- 
lution equations of the system. 



The present paper is organized as follows. The ba- 
sic equations describing flnite temperature Bose-Einstein 
condensates, as well as the properties of the condensate 
dark matter particles, are reviewed in Section |TT1 The 
equations of state of the finite temperature dark mat- 
ter, trapped by a gravitational potential, are obtained in 
Section Hm The cosmological evolution of the finite tem- 
perature condensed dark matter is considered in Section 
rvl We discuss and conclude our results in Section FVl 



II. FINITE TEMPERATURE BOSE- EINSTEIN 
CONDENSATE DARK MATTER 

In the present Section we briefly review the basic prop- 
erties of the high temperature Bose-Einstein conden- 
sate, and of the dark matter condensed particles. For 
a detailed discussion of the considered issues and of the 
derivation of the main results we refer the reader to the 
papers and books in [23[-|27J. 



A. The generalized Gross-Pitaevskii equation and 
the hydrodynamic representation 

The Hciscnbcrg equation of motion for the quantum 
field operator $ describing the dynamics of a Bose- 
Einstein condensate at arbitrary temperatures is given 
by mH 



ih 



9$ (r, t) 



2m 



A + mVgrav (r, t) + 



.g'<l+ (r, t) $ (f, t) 



$(r,t), (1) 



where m is the mass of the condensed particle, Vgrav (r", t) 
is the gravitational trapping potential, and g' ~ 
ATTlafi^/m, with la the s-wave scattering length. Eq. dU) 
is obtained under the assumption that the interaction 
potential can be represented as a zero-range pseudo- 
potential of strength g' . By taking the average of Eq. ([T]) 
with respect to a broken symmetry non-equilibrium en- 
semble, in which the quantum field operator takes a non- 
zero expectation value, we obtain the evolution equation 

for the condensate wave- function ^ {f,t) = ( $ (r , t) \ . 

Hence for the exact equation of motion of ^ (?% t) we find 



ih 



a* (?% t) 
'di 



-7^— A + mVgrav {r, t) 
Zm 



^ (f, t) + 
g'($+(r,t)$(f,t)<|.(f,t)\. (2) 



By introducing the non-condensate field operator 
'4'{f,t) so that l>(f, t) — ^(r, f) -I- ■ip{r,t), where the 

average value of ij} (r, t) is zero, (ip{r,t)\ = 0, wc can 

separate out the condensate component of the quantum 
field operator to obtain the equation of motion for ^ as 

follows lilKii 



ih 



(9* (f, t) 
dt 



h^ 

- -z—A + mVgrav (r, t) + gpc (r, t) 
2m 



^9P{r,t) 



* (r, t) + gpm (r, t) ** (r, t) + 



9Pi,+i,i, (r, t) 



(3) 



where we have denoted 



9 



inlah^ 



m^ 



(4) 



and we have introduced the local condensate mass den- 
sity 

p^{r,t) = mnc{f,t) ^ m\^ ir,t)\\ (5) 

the non-condcnsatc mass density 

p (r, t) = mn (r , t) = m ( -0"*" (r, t) ip (r, t) ) , (6) 
the off-diagonal (anomalous) mass density 

Prn {r, t) = mm (f, t) = m N) {r,t) 'ij) {r,t)\ , (7) 
and the three-field correlation function density 

respectively. 

In the following wc will restrict our analysis to the 
range of finite temperatures where the dominant ther- 
mal excitations can be approximated as high energy non- 
condensed particles moving in a self-consistent Hartree- 
Fock mean field, with local energy J22|425J 

£p {r, *) = ^ + mVgrav (r, t) + 2g [p^ (f, t) + p (f, t)] = 



f 



j^ + U.ffir,t), 



(9) 



where Ueff {r, t) = mVgrav {r, t) + 2g [pc (r, t) + J3 (r , t)] . 
Therefore in the present approximation we neglect the 
mean field effects associated with the anomalous den- 
sity prh and and with the three-field correlation function 

/■!/j+'0-0 y In the case of dark matter halos with a large 

number of particles this represents a very good approx- 
imation, the contribution of the anomalous density and 
of the three-field correlation function to the total density 
being of the order of a few percents [23] ■ 

In the thermal cloud the collision between particles 
forces a non-equilibrium distribution to evolve to the 
static Bose-Einstein distribution /" (r , p) [221 ■ Hence the 
particles in the thermal cloud are in thermodynamic equi- 
librium among themselves. By using a single-particle rep- 
resentation spectrum the equilibrium distribution of the 
thermal cloud can be written as 



fip,r,t) 



^Pep{r,t)-jl 



1 



(10) 



where /3 = l/fcsT, with ks Boltzmann's constant, and p 
is the chemical potential of the thermal cloud. In order 
to determine p we assume that the condensate and the 
thermal cloud components are in local diffusive equilib- 
rium with respect to each other. The requirement of a 



diffusive equilibrium between the cloud and the conden- 
sate imposes the condition |22| - [24j 



^J■c = M, 



fill 



where ^c is the chemical potential of the condensate. 
Therefore /ij, also determines the static equilibrium dis- 
tribution of the particles in the thermal cloud. 

The equilibrium density of the thermal excitations is 
obtained by integrating the equilibrium Bose - Einstein 
distribution over the momentum. Thus we obtain |22h25I 



(^^' *) = TTTTs / d'pf (P, r, t) = -^53/2 [z (f, t)] 
(zTrn.) J ■^T 



m 
T 



where At = ^2'Kh?l3/m is the de Broglie thermal wave- 
length, 53/2(2:) is a Bose-Einstein function, and the fu- 
gacity z (f, t) is defined as 

z (f, t) = e/3[A-c/.//(^-t)] = ^-Pgp.(P,t)_ (^3) 

The pressure p of the thermal excitations can be ob- 
tained from the definition u2, Sjl 



P (r, t) = 
and is given by [2^, l24 
P {r, t) = 



dp p^ 
{2Trhf 3m 



f{p,r,t).. 



-TT3-95/2 [z (r, t)] . 



(14) 



(15) 



The generalized Gross-Pitaevskii equation can be 
transformed to a hydrodynamic form by introducing 
the Madelung representation of the wave function as 
^ (i, r) = y^exp [{i/h) S (r, t)] . Then by taking into ac- 
count that in the present approach we neglect the effects 
of the mean field associated with the anomalous density 
and the three-field correlation function, it follows that 
Eq. ([3]) is equivalent to the following hydrodynamic type 
system |22h24| 



dpc 
dt 

dS__ 



V • {PcVc 



0, 



(16) 



(17) 



where the local velocity of the condensate is given by 
Vc {r, t) = {h/m) VS. The chemical potential of the con- 
densate is defined as 



Mc 



2m ,/rt7 



+ mVgrav {r, t) + gpc {r, t) + 2gp [r, t) . 

(18) 
Eq. (fT7|) can be reformulated as the Euler equation of 
fluid dynamics for the condensate. 



dvc 

i—r- 

dt 



dvc 
dt 



+ [Vc ■ V) Vc 



= -Vlic- 



(19) 



B. Physical properties of the dark matter 
condensate particle 

The zero-temperature approximation, which gives a 
very good description of the galactic dark matter halos 
[ijj 113j allows us to make an estimate of the physical 
properties of the dark matter particle. From the analysis 
of the static Bose-Einstein condensate dark matter halos 
it follows that the radius R of the co ndensate dark matter 
halo is given by i? = iT\Jh^la/Grn^ [IJ]. The total mass 
of the condensate dark matter halo M can be obtained 
as M = A-k'^ [h^la/Grrt') pc = AR^pgh/ir, where pgh 
is the central density of the galactic halo, giving for the 
mean value < p > of the condensate density the expres- 
sion < p >= Spgh/TT^- The dark matter particle mass in 
the condensate is given by [IJ] 



m 



~gW 



1/3 



6.73 X 10" 



1I/3 



\la{irn)Y"[R{k^c)]-^'^ eV. 



(20) 



For Za « 1 fm and _R w 10 kpc, the typical mass of the 
condensate particle is of the order of m k, 14 mcV. For 
la ~ 10® fni, corresponding to the values of la observed 
in terrestrial laboratory experiments, m sa 1.44 cV. 

An important method of observationally obtaining the 
properties of dark matter is the study of the collisions be- 
tween clusters of galaxies, like the bullet cluster (IE 0657- 
56) and the baby bullet (MACSJ0025-12). From these 
studies one can obtain constraints on the physical prop- 
erties of dark matter, such as its interaction cross-section 
with baryonic matter, and the dark matter-dark matter 
self-interaction cross section. If the ratio a^n = ct /m of 
the self-interaction cross section a = Awla and of the dark 
matter particle mass m is known from observations, then 
the mass of the dark matter particle in the Bose-Einstein 
condensate can be obtained from Eq. (f20|) as [23| 



^3/2 ft2 



2G i?2 



2/5 



(21) 



By comparing results from X-ray, strong lensing, weak 
lensing, and optical observations with numerical simula- 
tions of the merging galaxy cluster IE 0657-56 (the Bul- 
let cluster), an upper limit (68 % confidence) for am of 
the order of am < 1-25 cm^/g was obtained in [28|. By 
adopting for am a value of am = 1.25 cm^/g, we obtain 
for the mass of the dark matter particle an upper limit 
of the order 



TO < 3.1933 X 10^3^ 



R 



10 kpc 



-4/5 



1.25 cm2/g 



1.25 cm2/g 



1/5 



g = 0.1791 X 



R 



10 kpc 



-4/5 



1/5 



meV. 



(22) 



By using this value of the particle mass we can estimate 
the scattering length la as 



density of the thermal cloud can be represented as 



la < 



Om X m 



47r 



1.7827X 10"^^ cm = 1.7827 x 10"^ fm. 



(23) 

This value of the scattering length la , obtained from the 
observations of the Bullet cluster IE 0657-56, is much 
smaller than the value of Z^ = 10'* — 10^ fm correspond- 
ing to the EEC's obtained in laboratory terrestrial ex- 
periments [lll |. 

A stronger constraint for Om, was proposed in [29[, so 
that Um G (0.00335 cm^/g, 0.0559 cm^/g), giving a dark 
matter particle mass of the order 



(9.516 X 10"^^ - 1.670 X 10"^^) 



R 



lOkpc 



-4/5 



g = 



(0.053-0.093) 



R 



lOkpc 



-4/5 



meV, 



and a scattering length of the order of 



(5.038 
(5.038 - 



27.255) X 10"^^ 
27.255) X 10 



cm 
fm. 



(25) 



III. 



EQUATION OF STATE OF THE FINITE 
TEMPERATURE BOSE-EINSTEIN 
CONDENSATE 



In the present Section we obtain the equation of state 
of the Bose-Einstein condensate dark matter in thermal 
equilibrium with the thermal excitations described as a 
non-condensate cloud. To obtain the equation of state 
of the global condensate plus thermal excitations cloud 
we proceed in two steps: first we explicitly construct the 
equation of state of the thermal excitations, and then 
we obtain the temperature dependent equation of state 
of the condensate. The thermodynamic parameters of 
the condensate plus thermal excitations system can be 
obtained by adding the densities and pressures of the 
thermal excitations and of the condensate, respectively. 



A. The equation of state of the thermal excitations 

In order to obtain an easy to handle form of the matter 
density of the finite temperature non-condensate particle 
in thermodynamic equilibrium with the condensate, we 
power expand the .93/2 (^) Bosc-Einstcin function, so that 

c,3/2(e-^) = 2.612 - 3. 544V^ +1.4602;- 

0.103a;2 -^ 0.00424a;3 -F O fx^/^j . (26) 

For a; < 1 Eq. (I^5|) approximates the function 
53/2 (g"^) with an error smaller than 1%. Therefore the 



P = 



\rp 



2.612 - 3.544 V^gp^ + 1.460/3.g/9c - 



0.103 (/35p,)'+ 0.00424 (/3gp,)'' 



(27) 



The previous results can be written in a more transpar- 
ent form if we introduce the condensation temperature 
Ttr. given by [H 



Tt, 



o fc2 2/3 



C2/3(3/2)m5/3fcB' 



(28) 



(24) 



where ^(3/2) is the Riemann zeta function, and pt,. is the 
density of the dark matter at the condensation moment. 
Then we obtain immediately 



m 



Ptr 



C(3/2) 



T 

Ttr 



3/2 



and 



/3,gp, = 2(2/3(3/2) 



T 



1/3 2/3 \ T, 



Pc, 



(29) 



(30) 



respectively. By introducing the dimensionless conden- 
sate density 9, related to the condensate density pc by 
the relation 



Pt 



2/3"! 



1/3 



la ' 

we obtain the thermal cloud density as 



(31) 



PiT,0) = Ptr 



T 

Ttr 



3/2 



1 - 2.642 



T 

Ttr 



-1/2 



V~e 



T 
2.120 ( — 

. Ttr 



T 
0.572 ( —- 

. Ttr 



0.088 



T 

Ttr 



(32) 



At the initial moment of the condensation p^ = 0, and 
therefore 9 {Ttr) = 0. An order of magnitude estimate 
of the maximum value of 9 can be obtained by assuming 
that the density of the condensate is of the same order of 
magnitude as the initial transition density, pc « ptr, thus 
giving 0„,ax ~ lap\i^/m^'^. For Ttr - 10*° K, we obtain 
Ptr = 7.262 X 10-22 g/cm^, and ^^ax = 1-936 x lO'^^. 
However, due to the incertitude in the numerical values 
of the physical parameters, we will consider the behavior 
of the condensed dark matter system for larger values of 
^max- The variation of the thermal cloud density p is 
represented, as a function of 6, and for different values of 
T/Tt„ in Fig. m 
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FIG. 1: The density of the thermal cloud p as a function of 
e for different values of T /Ttr- T/Ttr = 0.95 (solid curve), 
T/Ttr = 0.90 (dotted curve), T/Ttr = 0.85 (dashed curve) 
and T/Ttr = 0.80 (long dashed curve), respectively. 



FIG. 2: The pressure p/po of the thermal cloud as a function 
of e for different values of T/Ttr- T/Tr = 0.95 (solid curve), 
T/Ttr = 0.90 (dotted curve), T/Tr = 0.85 (dashed curve) 
and T/Ttr ~ 0.80 (long dashed curve), respectively. 



To obtain the pressure of the non-condensate particles 
we start with the expansion 

£,5/2 (e-^) = L341 + 2.363x3/2 -2.612X- 

0.730a;2 + O.OMGx^ + O (x^^^) , (33) 



valid for a; < 1, which gives 



p 



1.341 + 2.363 {Pgpc) ' ~ 2M2Pgpc 



3/2 



Arp 



0.730 iPgp.f + 0.0346 {pgpcf 



In terms of the transition temperature we have 



ksT _ ksTtr 

-IT- - Ptr 

Arp 



T 



C (3/2) 771 \Ttr 



5/2 



(34) 



(35) 



and thus we obtain for the non-condensate pressure the 
expression 



p{T,e) = pt, 



ksTtr 



T 

Ttr 



0.513 4-6.684 



5/2 



T 

Ttr 



-3/2 



g3/2 



3.793 ( — ) 6*- 4.022 ( — ' ^^ 

, Ttr J V '^tr , 



0.724 



T 

Ttr 



-3 



(36) 



The variation of the thermal cloud pressure p/po; 
where po = PtrkBTtr/'m, is represented, as a function 
of 6*, in Fig. [1 

The equation of state p = p{p) of the thermal cloud is 
represented in Fig. [31 
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FIG. 3: The pressure p/po of the thermal cloud as a function 
of p/ptr for different values of T/Ttr'. T/Tr = 0.95 (solid 
curve), T/Ttr = 0.90 (dotted curve), T/Tr = 0.85 (dashed 
curve) and T/Tr = 0.80 (long dashed curve), respectively. 



The obtained equations of state of the thermal excita- 
tions arc valid only in the temperature range determined 
by the condition 



or, equivalently, 
T 

Ttr 



ksT > gpc 



> 2(2/3(3/2) ^ 



1/3 2/3' 



(37) 



(38) 



The equation of state of the finite temperature 
condensate 



In order to obtain the equation of state of the con- 
densate we introduce the Thomas-Fermi a ppr oximation 
for the condensate wave function |22| - [23. In the 



Thomas-Fermi approximation, the kinetie energy term 
— {h^ /2m) A of the condensate particles is neglected. 
Hence in this approximation the chemical potential of 
the condensate is given by [2^ 



Mc = mVgrav (r, t) + gpc {r, t) + 2gp (f, t) , 



(39) 



Eq. PT)) can be reformulated as the Euler equation of 
fluid dynamics for the condensate, 



Pc 



dvc 
'dt 



Pc'^Vgrav PcV [pc + 2p\ 



(40) 



By taking into account the identity p^ p'l = 
V [n/ (n + 1)] pJ?+\ and with the use of Eq. d^, it fol- 
lows that Eq. (PO)) of the motion of the Bosc-Einstcin 
condensate can be written as 



dvc 
'~dt 



-Pc^V, 



grav 



"^Pc 



(41) 



where the pressure of the finite temperature Bosc- 
Einstein condensate in thermal equilibrium with a gas 
of thermal excitations is given by 



Pc {T, Pc) 



5 2 „ 2.362^ (fc^r)-^/^^ + 

„2 3 

1.460f3- [keT)'^ pi - O.lSTfg- x 



Xrp 



\rp 



{kBT)-'pl 



0.00636^3- (fcsT) 



3 „4 



A3 --^ ' P' 



(42) 



In terms of the dimensionless density 0, defined by 
Eq. ((3T|) . we obtain 



Pc{T,9) = pt, 
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FIG. 4: The normalized pressure Pc/po of the condensate as a 
function of for k = 50 and for different values of T/Ttr'- 
T/Ttr = 0.95 (solid curve), T/Ttr = 0.70 (dotted curve), 
T/Ttr = 0.60 (dashed curve) and T /Ttr = 0.40 (long dashed 
curve), respectively. 



metric, given by 

ds^ = -c^dt^ + a^{t) {dx^ + dy^ + dz^) , (45) 

where a{t) is the scale factor describing the cosmological 
expansion. The Hubble function H{t) is defined as H = 
a/a. 

For the matter energy-momentum tensor we consider 
the case of the perfect fluid energy-momentum tensor, 
given by 



T'^" = {ptotc'^Ptot)u^vy +Ptotg^\ 



(46) 



(43) 



The ratio Pc/po for the condensate is determined by the 
numerical value of the dimensionless parameter 



where ptot and ptot represent the total energy density and 
pressure of the Universe. 

As for the matter content of the Universe, we assume 
that it consists of radiation, with energy density prad 
and pressure Prad^ prcssurcless {pb = 0) baryonic matter, 
with energy density pb, dark matter, with energy density 
Pdm, and pressure pdm, respectively, and dark energy, 
described by the cosmological constant A. Hence the 
total energy density and pressure of the matter in the 
Universe is given by 



^1/3 

1/3, • 
Ptr 'a 



(44) 



and 



The variation of the condensate pressure with respect 
to 9 is represented, for k, = 50, Fig. |4l 



IV. COSMOLOGICAL EVOLUTION OF THE 

UNIVERSE WITH FINITE TEMPERATURE BEC 

DARK MATTER 

We assume that the space-time geometry of the Uni- 
verse is the flat Friedmann-Robcrtson- Walker (FRW) 



Ptot = Prad + PDt,[ + Pb, 



Ptot = Prad +PDM, 



(47) 



(48) 



respectively. In the following we neglect any possible 
interaction between these components, by assuming that 
the energy of each component is individually conserved. 
Thus, the gravitational field equations, corresponding to 
the line element given by Eq. (j45p . become 



3— = SttG (pb -I- Prad + Pdm) + A, 
a'' 



(49) 



BttG 



(Prad +PDm) + a, 



(50) 



P^ + i (p^ + ^)- = 0, i^b, rad, DM. (51) 



From Eq. ((5T|) it follows that the energy density of the 
radiation and of the baryonic matter are given by 



and 



Prad 



Pb 



Prad,0 



(a/flo) 



Pbfl 

{a/aaf 



4' 



(52) 



(53) 



respectively, where Prad,o and pt^ are the mass densities 
corresponding to a = cq. 



A. Evolution equations of the tw^o-component finite 
temperature BEC darlt matter 

Even if the total energy-density of the condensate- 
thermal excitations system is conserved, the energies of 
the condensate and of the thermal cloud arc not inde- 
pendently conserved, due to the growth of the conden- 
sate through absorption of the particles from the thermal 
cloud. This process must naturally occur during the ex- 
pansion and subsequent cooling of the Universe. More- 
over, the collisions in the thermal cloud will transfer a 
particle to the condensate, so that the number N of the 
condensed particles will increase as N — >■ N + 1, and 
there is of course the reverse process, where a collision of 
a thermal cloud particle with one within the condensate 
transfers a particle from the condensate into the ther- 
mal excitation band, so that N ^ N — 1. However, in 
the following we will consider only a simplified model, 
in which the particles of the thermal cloud turn into the 
condensate at a rate 1/tcoI particles per second, lead- 
ing to a decrease of a number of particles in the thermal 
cloud and to the corresponding increase of the number of 
the condensate particles. Therefore the equations of the 
energy conservation of the two dark matter components 
take the form 



— + 3- [~+-^ 
dt a \ c^ 



and 



dpc 
dt 



a V 



= -r™*(Ar), 



r"(iv). 



(54) 



(55) 



respectively, where r°"*(iV) and F'" describe the time 
variation of the density of the two components in the 
finite temperature dark matter system. From Eqs. (j54p 
and (|55p it follows that the total dark matter density 
PDA! = P + Pc satisfies the conservation equation 



dpDM 

dt 



3-(, 



Pdm 



PPM \ 



0, 



(56) 



where pdm ~ P + Pc is the total pressure of the 
condensate-thermal excitations dark matter system. 

The time rate TcoI at which the excited components 
turn into the condensate may be approximated by the 
transition probability W{N) that can be estimated by 
using the quantum kinetic theory, so that 1/tcoI ~ W{N) 
[SO]. W{N) gives the rate at which the thermal particles 
above the condensate energy band enter the condensate 
due to inter-particle collisions. The main assumption in 
obtaining W{N) is that the condensate does not readily 
act back on the thermal component to change its tem- 
perature. This assumption is expected to be valid when 
one considers equilibrium or quasi-equilibrium situations. 
The transition rate of the thermal excitations to the con- 
densed state can be given as [SO] 



Tool T^K^ 



Pc 



knT 



Ki 



Pc 



knT 



(57) 



where jl and /i^ arc the chemical potentials of the thermal 
cloud and of the condensate, and A'l is a modified Bessel 
function. Eq. ([57]) can be simplified to the form 



1 



a, 



4m {IgkBTf 



WirhCr 



T 



Tcoi " Tin- C^/3(3/2)m?2 \Ttr, 

(58) 
where the correction factor Ccon assumed to be a con- 
stant, incorporates the effects of the exponential and of 
the Bessel function, respectively. The transition rate can 
be written as 



1 

Tcol 



Crnr X 2.0705 X 10 



'28 



10 



-20 



cm 



10-37 g^ 

(59) 



Generally, the source terms in the thermal cloud and 
continuity equations can be written as [23 



crpc 



■i ,r. ^73/ dp2dpsdpi5 {pc +P2-P3- Pi) X 

TTTTi'' (ZTTh) J 

S (ec + £2 - £3 - £4) X /2 (1 + h) (1 + U) , (60) 



and 



crp 



- — -3 / dp2dp3dpi5 {pc +P3-P2- Pi) X 
■nh) J 



TTm^ (2'Kh) 

'5(£c + £3-£2-e4)x/2(l + /3)/4, (61) 

where u is the scattering cross section, p and e denote the 
particle momentum and energy, and fi, i = 1, 2,3,4 are 
the Bose-Einstcin distribution functions of the thermal 
particles and of the condensate particles, respectively. 
The out collision rate represents the scattering of an in- 
coming thermal particle (2) from the condensate to pro- 
duce two thermal atoms (3) and (4). The designation 
"out" means that a particle is leaving the condensate. 
Thus, r°"* is the rate of increase of the number of ther- 
mal particles per unit volume and per unit time as a 



result of the collision with a condensate particle. The 
reverse process gives the " in" collision rate of a thermal 
particle F™. In the thermodynamic equilibrium state 



r, 



(62) 



and in the following we assume that this condition holds 
for all times. The equality between the transition rates 
between the thermal clod and the condensate also allows 
us to define the equilibrium collision time t^oi, given by 
Eq. dSHl). Generally, r°"* and T*" can be obtained by 
numerical integration of Eqs. ([501 and (pT|) . However, 
in order to obtain a semi-analytical model of the cosmo- 
logical evolution of the finite temperature Bosc-Einstcin 
condensate dark matter, for the energy density trans- 
fer rate T we introduce the simplifying phenomenological 
assumption that it is proportional to the density of the 
condensate pc, and we write it as 

T[N) ex pc = AK — pc = 

Tcoi T^ri^ 



WnhK 



T 



C4/3(3/2)mZ2 \Ttr 



1.472 X lO'^^ X K X 



Pc = 



10-2" cm 



10-37 g 

2 



T 



,(63) 



where K is a constant. In terms of the dimensionless 
condensate density we obtain 

\mrni... „ / / \ 

r 



C4/3(3/2)m2 

WWm 



Kpl ( — 



T 

Ttr 

T 

Ttr 



(64) 



B. 



Cosmological dynamics of the finite 
temperature BEC dark matter 



In the following we denote the " reduced" dimensionless 
temperature T/Ttr by Xj X = T/Ttr- By taking into 
account that in the case of the finite temperature Bose- 
Einstein condensate dark matter the energy densities pc 
and p are functions of the reduced temperature x ^-nd of 
the condensate density 9, pc = pc (x, ^), P ^ P (Xi ^)i the 
energy conservation equation for the dark matter takes 
the form 



dep{x,0)d + d^p{x,e)x + 



p{x,e) 



and 



3^ 
a 



p{x,o) 



{x,0)9 + d^Pc{x,0)x 
Pc{x,0) 



r. 



(65) 



Pcix,9) 



= r, 



(66) 



where de = d/dO and d^ = d/dx, respectively. 

In the present Subsection we consider the cosmological 
evolution of the condensate dark matter component only, 
by ignoring the effects of the radiation and of the baryonic 
matter. In order to simplify the formalism we introduce 
the critical density of the Universe as pcr = Hq/8ttG, 
where Hq is the value of the Hubble function at a ~ uq, 
and we define the density parameter of the dark matter 
as ^DMdr ~ Ptr/ Per- We also rcscale the comoving time 
variable t by introducing the dimensionless time variable 



''' — HQy'Q£iM,trt, 



(67) 



respectively. Accordingly, the constant K in the expres- 
sion of r is rescaled as ii" = K/Ho^/TljxM^, and we 
denote the corresponding energy density transition rate 
by r. Moreover, we normalize the scale factor so that 
ao — 1 gives the value of the scale factor corresponding 
to the present age of the Universe. 

By solving Eqs. ([HS]) and (|66l) for x and 9, we obtain 
the following differential equations describing the time 
variation of the temperature and of the condensate frac- 
tion. 



dx 
dT 



and 



i \ da 
A a dr 



P 



depc - p. 



Pc 



dep 



A 



{dePc + dep) . 



(68) 



d0_ 



A a dr 



(/'<= + 5) ^xP + ^^xP' (69) 



respectively, where we have denoted A = depcd^p — 
dxPcdep = dopcd^p, and we have taken into account 
that Pc is independent of the reduced temperature Xi 
dxPc = 0. 

The evolution of the scale factor of the Universe is 
described by the equation 



1 da 
a dr 



I Pdm{x,&) 



Ptr 



(70) 



Eqs. ([68|) . (|69|) and ((70)) give a complete description of 
the cosmological evolution of the finite temperature Bose- 
Einstein condensed dark matter. In order to completely 
determine the cosmological dynamics of the dark matter, 
we need the initial value of 9 and T at a = atr, 9 (atr) = 
9o. Since at the initial moment of the transition there 
is no dark matter in the form of a condensate, it follows 
that 6*0 = at a = Otr- As for the initial value of the 
reduced temperature, it is given by xio'tr) = 1- In the 
following we denote 



Ko 






(71) 



The cosmological evolution of the Bosc-Einstein con- 
densed dark matter depends on two dimensionless pa- 
rameters, kq and K. By assuming for the condensation 



10 



temperature Tjy. a value of the order of Tfr = 10 K reaches a maximum value ^max at r 



But for 



[20|, for a particle mass of the order of m = 10^''' g, we 
obtain kq ~ 1.533 x 10^°. The transition density cor- 
responding to this temperature is ptr = 7.255 x 10^^^ 
g/cm^, giving, for la = lO^^o ^^^^ ^ ^ 5^53 ^ ioi4. Por 
a transition density of the order of ptr = 10^^ g/cm^, 
and for la ~ 10~^° cm, we obtain kq = 1.897 x 10^^ and 
K — 4.641 X 10^". However, due to the poor knowledge 
of the transition temperature and density, the values of 
these parameters are very uncertain. That's why in the 
following we will consider a wide range of values for the 
parameters kq and k, ranging from small to high values. 
The constants kq and k are not independent, but they 
are related by the relation 

27rft2 



1 1 



C2/3(3/2)c2 7712^2^2 



4.088 X 10 



39 



10 



-37 



la 



10 



-20 



1 



(72) 



Once the values of k and kq arc fixed, or obtained, for 
example, from general physical/observational considera- 
tions, the value of the product mla is determined as 



mla = 



/2^ 



1 



6.394 X 10 



-38 



Ci/3(3/2)cKy^ 
In the following we also denote 



KwKq 



g cm. (73) 



K\ = K/la. (74) 

The time variation of the dimensionless condensate 
density is represented in Fig. \5\ 



o 




8.x 10" 



1.x 10" 



time intervals so that r > Tmax, the cosmological ex- 
pansion takes over, and the density of the Bose-Einstein 
condensate dark matter decreases in time. The variation 
of the temperature of the condensate dark matter after 
the start of the phase transition is represented in Fig. [51 




FIG. 6: Time variation of the reduced temperature x = T/Ttr 
of the finite temperature condensed dark matter, for Ki = 5 x 
10^^ (solid curve), Ki = 4x 10^'' (dotted curve), Ki = 2x 10^^ 
(dot dashed curve), and Ki = 0.6 x 10'^^ (dashed curve), 
respectively. The numerical values of the initial conditions 
and of the parameters are 9q = 10~^, atr = 10~*, xo = 
0.9999, K. = 10^^ and kq = 10^". 



The time variation of the total density of the dark mat- 
ter, as well as of its total pressure, are represented in 
Figs. [7] and [51 respectively. 



s 

□ 

o 

O 
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2.X10"" ^xlO"'"* e.xlO"'"* S.xlO"" l.xlO"'^ 



FIG. 5: Time variation of the dimensionless condensate den- 
sity e for A'l = 5 X lO^'^ (solid curve), K^^Ax 10^^ (dotted 
curve), JsTi = 2 X 10^^ (dot dashed curve), and Ki = 0.6 x 10^'' 
(dashed curve), respectively. The numerical values of the 
initial conditions and of the parameters are 9o = 10~^, 
atr = 10-^, xo = 0.9999, k = 10^' and kq = 10^°. 

Due to the decrease of the temperature during the cos- 
mological evolution, the condensate density increases im- 
mediately after the condensation phase transition, and it 



FIG. 7: Time variation of the total density poM of the fi- 
nite temperature condensed dark matter, for Ki = 5 x 10^^ 
(solid curve), _?Ci = 4 x 10^'^ (dotted curve), A"i = 2 x 10^'' 
(dot dashed curve), and Ki — 0.6 x 10"^^ (dashed curve), re- 
spectively. The numerical values of the initial conditions and 
of the parameters are 60 = 10~^, atr ~ 10"*, xo ~ 0.9999, 
K — 10^^ and Ko = 10^". po is the value of the condensate 
density corresponding to 6q and xo. 
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FIG. 8: Time variation of the total pressure pdm of the fi- 
nite temperature condensed dark matter, for Ki = 5 x 10''^ 
(solid curve), iCi = 4 x 10^'' (dotted curve), Ki ^ 2 x 10^'' 
(dot dashed curve), and Ki = 0.6 x 10'^^ (dashed curve), re- 
spectively. The numerical values of the initial conditions and 
of the parameters are 9o — 10~^, atr = 10"'*, xo ~ 0.9999, 
K = 10*^ and Ko = 10*". 



described by the scale factor a^"''^\ From Eqs. (|49)) and 
(|5ip we obtain immediately 



Pdm _ fatr^^ 

Ptr 



(v)' 



and 



2/3 



«^"^nT)-U \r + ^<'r. 



3/2 



2/3 



(75) 



(76) 



respectively. The comparison of the time variation of the 
scale factor of the Universe filled with condensed dark 
matter with the standard pressureless ACDM dark mat- 
ter model is represented in Fig. [S] 

As one can see from the Figure, in the presence of the 
condensed dark matter the expansion of the Universe is 
faster. 



C. Cosmological evolution of the Universe with 
BEC dark matter, dark energy and radiation 



We compare the cosmological evolution of the finite 
temperature condensed dark matter with the standard 
ACDM model of the pressureless dark matter, with to- 
tal density p^^j, and with the cosmological expansion 



In the case of a Universe filled with dark energy, radia- 
tion, baryonic matter with negligible pressure, and Bose- 
Einstein condensed dark matter, respectively, the time 
evolution of the scale factor is given by the differential 
equation 



1 da 
a dt 



H{)\ 



In 



b.O 






id,0 



n 



DM,tr- 



Pdm 

Ptr 



VL\A > tt, 



(77) 



which must be integrated with the initial condition 
o-{ttr) = oitr- In the following for the Hubble constant 
we adopt the value H^ = 70 km/s/Mpc = 2.273 x 10"*^ 
s~^ |3l| . giving for the critical density a value of Pco = 
9.248 X 10~3" g/cm^. For VtoMfi ^b,o, ^rad,o, and Vt^ 
we adopt the numerical values r^DM.o = 0.228, iTh^o = 
0.0456, nradfi = 8.24 x lO'^, and 17a = 0.726 [sj', re- 
spectively. By taking into account that ptrO-tr ~ Pc,q, 
it follows that 0.£iM,tr ~ ^_DM,o/«tr) ^nd the evolution 
equation of the scale factor can be written as 



1 da 
a dr 






1 



' »fc,0 1 

^DALtr a'^ nDM.tr a* 



Pdm 

Ptr 



n^ 



n 



DM.tr 



(78) 
The comparison of the time variation of the scale factor 
of the Universe filled with dark energy, radiation, bary- 
onic matter with negligible pressure, and Bose-Einstein 
condensed dark matter, respectively, and with the stan- 
dard pressureless ACDM dark matter is represented in 
Fig.Hni 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the cosmolog- 
ical dynamics of the finite temperature Bosc - Einstein 
Condensed dark matter. In a simplified approach one 
can describe a finite temperature condensed system as a 
mixture of two fluids in thermodynamic equilibrium, the 
thermal cloud, consisting of the thermal excitations of the 
condensate, and of the proper condensate state. At the 
beginning of the condensation process the bosonic dark 
matter is mostly in the form of the thermal cloud. How- 
ever, once the temperature of the system drops, more and 
more particles from the thermal cloud will enter into the 
condensate energy band. Therefore, the cooling of the 
Universe due to the cosmological expansion did lead to 
the significant decrease in the number of particles in the 
thermal cloud, and to the corresponding increase in the 
number of the condensate particles. We have described 
the cosmological transition process in the framework of 
a two interacting fiuid model. Due to the cooling of the 
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FIG. 9: Time variation of the scale factor of the finite tem- 



perature condensed dark matter for K\ = 5 x 10' 



,37 



(dotted curve), Ki = 2 x W 

l37 



(solid 
(short- 



curve), /sTi = 4 X 10^^ 

dashed curve), K\ = 0.6 x 10^^ (long-dashed curve), and of 
the pressureless dark matter ACDM model (dot-dashed). The 
numerical values of the initial conditions and of the physical 
parameters are Oo = 10"^, atr = 10"", xo = 0.9999, k = 10^^ 
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pens at the moment when the time rate of the decrease 
in the condensate density due to the expansion of the 
Universe becomes higher than the density transfer rate 
from the thermal cloud. As one can see from Eq. ([M)) 
the maximum value of the condensate density is reached 
■when dO/dr ~ 0, ■which gives the algebraic equation 



Ptr 



(79) 



a function of the reduced temperature x and of the physi- 
cal parameters of the condensed dark matter system. The 
contour plot of Eq. ([7^ is represented in Fig. [TTJ A con- 
tour plot is a plot of equipotential curves z{x,y), where 
each of the curve is the geometric locus of pairs {x, y} 
such that z = constant. In the case discussed here, each 
curve tracks the set of pairs {Omax, x} for which the value 
of the right hand side of Eq. ([7^ is a constant (the values 
of these constants are shown explicitly on each curve). 
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FIG. 10: Time variation of the scale factor of the Universe 
filled with dark energy, radiation, baryonic matter with neg- 
ligible pressure, and Bose-Einstein condensed dark matter 
for Ai = 4 X 10^'^ (solid curve), Ai = 3 x 10^^ (dotted 



^36 



-i36 



curve), Ai = 2 x 10'"' (short-dashed curve), Ai = 0.6 x 10'' 
(long-dashed curve) , and of the standard ACDM model (dot- 
dashed curve) . The numerical values of the initial conditions 
and of the physical parameters are Oq — 10~®, atr ~ 10"'', 
XO = 0.9999, K = 10^^ and kq = 10^". 



Universe, the thermal cloud looses its particles and the 
thermal excitations enter into the condensate. Therefore, 
the dynamics of the dark matter is determined by the 
increase of the density pc of the condensed component, 
which is initially increasing in time, and the correspond- 
ing decrease in the number of the thermal particles. How- 
ever, after reaching a maximum value, the density of the 
condensate component will start to decrease. This hap- 



FIG. 11: The contour plot of Eq. 
K= W^ andKo = 10"'. 



As a general result we have also found that the pres- 
ence of the condensed dark matter will accelerate the 
expansion of the Universe. In the case of zero temper- 
ature condensed dark matter this behavior has been al- 
ready pointed out |18l - [2(J |. However, a further increase 
in the speed of the expansion may be expected due to 
the finite temperature effects during the phase transi- 
tion. Moreover, due to the presence of the two interact- 
ing fluid components, the cosmological dynamics of the 
finite temperature condensed dark matter is much more 
complicated than expected. Therefore the two-fluid finite 
temperature condensation transition may provide some 
clear cosmological signatures that could help in discrim- 
inating between Bose - Einstein Condensed dark matter 
models, and the standard pressureless ACDM cosmolog- 
ical model. 
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